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We present in this paper transit-time variations in line-of-sight propaga- 
tion paths and systems operating at frequencies up to SO GHz. We discuss 
variations due to both atmospheric changes (no precipitation) and rain 
and point out some relationships to PCM systems. 

I. INTRODUCTION 

In n recent paper 1 , J. R. Pierce considered stable synchronization of 
large digital transmission networks, and pointed out that the realization 
of such a synchronized network calls for, among other things, more 
information concerning network transit-time variations. In this paper, 
we seek extreme values for these variations in line-of-sight propagation 
paths in order to provide some of this information. Estimates are given 
for (i) the maximum variation in transit-time, Ar mRX , which one might 
expect over the period of a year, and (u) the maximum time derivative, 
T mnx , which one might encounter. The estimated values are related to 
digital systems, with most specific examples given for a 500 megabit 
transmission rate. Variations due to changes in the atmosphere (no 
precipitation), and those due to rain are discussed separately. In the 
longer line-of-sight paths achieved in tandem systems, repeaters are 
assumed to be stable, that is, the concern herein is with atmospheric 
variations only. Delays associated with selective fading are not dis- 
cussed, but they are believed not to exceed the given estimates of the 
maximum variations. 

The transit-time t is given by the familiar relationship 



1 r p - 

= - nds (1) 

c Jp. 



where c is the velocity of light, n the medium refractive index, ds the 
differential path length, and the limits P, and P 2 represent the end points 
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of the path. If one assumes An max , the maximum index of refraction 
change which is expected over a given period of time, and if one assumes 
that the changes in n over the entire path S are perfectly correlated and 
equal to An max , then the corresponding maximum transit-time variation, 
AT msr , is given by 



At„ 



= - (A7W)£. 
c 



(2) 



In a subsequent section, AT m « will be computed after estimating An max . 
A similar set of assumptions concerning the time derivative n leads to 
an estimate of f max , namely, 



1 . 
= c W " 



t S. 



(3) 



II. MAGNITUDE OF TRANSIT-TIME VARIATIONS IN THE ATMOSPHERE 

Letting the superscript a denote the atmosphere, the estimate of 
AtC„ to be used in this section is based on data found in Bean and 
Dutton. 2 The data comprise eight years of point observations at six 
locations in the United States. The locations are listed in Table I where 
the maximum, minimum, and range An of the index of refraction are 
given in N units, with N = (n — 1) X 10 6 . A value which approximately 
strikes an average for the six locations listed will be used, namely, 
A?C X = 1-0 X 10~ 4 . Substituting in equation (2), 

ArLx = 0.333 X 10~ g S (4) 

where S is in kilometers. Equation (4) is plotted in Fig. 1. Since the 
atmospheric index of refraction, n a , may be regarded as independent of 
frequency up to 30 GHz 2 , the curve in Fig. 1 may be considered applic- 
able up to this frequency. 



Table I — Extremes in Index of Refraction for Six 
Locations in the United States 



Location 


max N 


min N 


AN 


Washington, D. C. 
San Antonio, Texas 
Bismarck, N. D. 
Colorado Springs, Col. 
Salt Lake City, Utah 
Tatoosh Island, Wash. 


393 
387 
368 
307 
323 
354 


277 
258 
260 
214 
230 
292 


116 

129 

108 

93 

93 

62 
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Fig. 1 — Estimated maximum atmospheric (no precipitation) transit-time varia- 
tions for long path length line-of-sight tropospheric communications systems. 

For long paths, At£, bx is overestimated because the assumption that 
An is perfectly correlated and equal over the entire path becomes un- 
realistic as S increases. However, the overestimation is not excessive 
because those components in An which primarily are due to seasonal 
and diurnal variations are highly correlated and combine to be of the 
order of 50 N units (approximately 30 N units for seasonal and 20 N 
units for diurnal) . Thus, Ar^ ax is conservative to within a factor of two 
for the continental United States. As an example, for a 3000-mile path 
in Fig. 1, At^ ox = 1.6 us; certainly, it would be no less than 0.8 us. In 
terms of a 500 Mb system, Ar a mai = 1.6 /xs is equivalent to an 800 bit 
variation. 

For short path lengths, more detailed data are needed and these are 
in Fig. 2. One notes that, for short paths, a An" which typifies a given 
region should be used instead of the A?i£, ax of 100 N units. Therefore, 
Fig. 2 shows AT m „ versus path length for each of the six locations listed 
in Table I. 

Now consider the problem of synchronizing two clocks separated by 
a distance S. For a digital system of pulse spacing T, two clocks may be 
considered synchronized if they are in phase to within a factor / of a 
pulse spacing. For the sake of argument in this discussion, we choose 
/ = 0.1, and we draw on Fig. 2 the lines /• T for pulse rates of 50, 100, 
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Fig. 2 — Maximum atmospheric (no precipitation) transit-time variations for 
short paths. Variations are shown for six locations in the United States. 

200, and 500 Mb. The intersections of the lines /• T with the curves for 
Ar mil « determine the maximum path lengths for which atmospheric 
transit-time variations may be neglected. For instance, these lengths 
are respectively 5, 2.5, 1.25 and 0.5 km for the 50, 100, 200, and 500 Mb 
rates when the region of concern is Washington, D. C. For some other 
criterion of synchronization, that is for some other /, these path lengths, 
of course, will be different. 

From Fig. 2, it is evident that for any but the shortest path lengths, 
some compensation for refractive index changes must be incorporated 
in a digital system running on a universal clock. Figures 1 and 2, how- 
ever, do not provide information on how rapidly one must compensate. 
This question is answered by first considering the time derivative of 
transit-time variations which are due to turbulence, and then by con- 
sidering those due to the motion of synoptic scale air masses. 



III. TIME DERIVATIVE OF TRANSIT-TIME VARIATIONS 

Table II shows the rms values of the hourly transit-time variations 3,4 
which are due to atmospheric turbulence. Because of the randomness 
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Table II — RMS Hourly Transit-Time Variations 
in the Turbulent Atmosphere 



RMS Delay 
(seconds) 


Path Length 
(miles) 


Frequency 
(GHz) 


1.5 X 10-' 2 

4 X 10~ 12 

15 X lO" 12 

22 X lO" 12 


2.25 
3.5 

10 

60 


10 
1 
1 
1 



of turbulent motion, one would expect the rms values to be proportional 
to the square root of the path length S. Thus, on an hourly basis, the 
variations for a transcontinental path work out to be 0.2 ns, which 
is one tenth of the pulse spacing in a 500 Mb system. On the basis of a 
five minute interval, measurements show that rms variations are two 
orders of magnitude smaller than the hourly changes. 4 Thus it appears 
that compensation for turbulence-induced atmospheric transit-time 
variations can be made on the order of tens of minutes for transcon- 
tinental links, and even more slowly for shorter paths. 

In contrast with the turbulence-related phenomena discussed above, 
the motion of synoptic-scale air masses brings about changes in index- 
of-refraction which are correlated over large regions, that is, for regions 
covering hundreds of kilometers, changes in n will be proportional to 
path length, S, not to (S)*. For instance, an advancing cold front brings 
with it a decrease in temperature 0, a decrease in water vapor pressure 
e, an increase in total pressure P, and an accompanying change AN F 
which is well correlated over the entire frontal advance. 

For the purpose of discussion, a model front is shown in Fig. 3. The 
model front passes transversely across a microwave transmission path. 
The index of refraction change AN F occurs over a distance S A , and the 




'-s A ->| 

Fig. '.i — Model of a frontal system moving across a line-of-sight transmission path. 
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front moves at velocity v. This model permits one to estimate the maxi- 
mum rate of change of transit-time f maI . 

Proceeding with the estimate of f max , from Fig. 3, 



= 10- AN " 



V = 



= 10 



_„A2V, 



<Si 'max '& 

where / A = S&/v. Substituting in equation (3), 

,-u AiV„ 



= 0.33 X 10" 



t- 



S. 



(5) 



Estimates of the quantity hN P /t L | max have been obtained from tem- 
perature, pressure, and relative humidity data continuously recorded 
at Crawford Hill since January 1, 1967. The records were examined for 
events of rapidly changing atmospheric conditions. N was computed 
for conditions just prior to and just after these events using the formula 2 



N = 77.6? + 3.73 X 10 5 4 



(6) 



where 6 is expressed in kelvins, and P and e in millibars. AN P was taken 
as the difference between the two computed values for each event. Table 
III gives the date, time, AN? , t& , and AN F /tt for the largest events 
recorded. It is noted from the Table that AN F /U |m» = -0788. Rounding 
off to .08, and substituting in equation (5), 



Table III — Index op Refraction Changes Accompanying 
Rapidly Varying Atmospheric Conditions 



r Date 


Time 


AN F 


t& (sec) 


ANp/tn 


2/28/67 


10:30 PM 


7.23 


900 


.00804 


7/14/67 


3:00 PM 


3.94 


600 


.00657 


10/3/67 


9:00 PM 


-13.67 


300 


.0456 


2/17/68 


2:00 PM 


13.63 


300 


.0455 


3/29/68 


3:30 PM 


15.93 


900 


.0177 


4/30/68 


6:00 PM 


7.55 


300 


.0252 


6/3/68 


4:00 PM 


13.75 


600 


.0229 


7/2/68 


8:00 PM 


22.18 


900 


.0246 


7/24/68 


2:30 PM 


19.84 


300 


.0662 


8/7/68 


2:30 AM 


18.67 


600 


.0311 


8/15/68 


4:00 AM 


-31.75 


1200 


.0264 


8/17/68 


3:30 AM 


6.18 


300 


.0206 


8/22/68 


9:40 PM 


-2.42 


300 


.00807 


11/29/68 


2:00 AM 


3.67 


180 


.0204 


12/5/68 


2:00 PM 


14.18 


180 


.0788 


6/13/69 


2:30 PM 


10.12 


600 


.0169 


6/24/69 


5:30 PM 


7.49 


300 


.0250 
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f max = 2.67 X 10" 13 5. (7) 

The maximum value S over which good correlation of N can occur will 
be taken as 500 km (this value will be discussed subsequently). From 
equation (7), T mnx ^ 1.36 X 10" 10 sec/sec. In terms of a 500 Mb system, 
T m . x = .0667 b/s, or l/16th of a bit per second. 

The derived T m „ is thought to be larger than would be encountered 
in practice. Complete parallelism (or coincidence), as shown in Fig. 
3, becomes highly improbable when the path length S approaches 500 
km. Primary reasons are: (i) a front and a transmission path can be 
oriented at angles ranging over a large fraction of ir radians; (ii) nature 
will not produce fronts composed of straight line segments, but rather 
of curves. Since a line-of-sight propagation route is made up of straight 
line segments, it will be improbable that a front would coincide with it. 
Consequently, it appears that the choice of S = 500 km is extreme, and 
that T max = 1/16 b/s at 500 Mb is a good upper bound. 

IV. TRANSIT-TIME VARIATIONS DUE TO RAIN 

Letting the superscript r denote rain, then An r is specified using com- 
putations of the medium refractive index for given rainfall rates under 
the assumption of a Laws & Parsons drop size distribution. 6 These com- 
putations include the index of refraction of the medium for 6, 16 and 
30 GHz for the rain fall rates 0.25, 1.25. 2.5, 5.0, 12.5, 25.0, 50.0, 100.0, 
and 150.0 mm/hr. For rainfall rates exceeding 150 mm/hr, proportional 
scaling of the 150 mm/hr medium index will be used. 

The nature of rainfall is such that higher fall rates are associated with 
smaller areas of coverage. For this reason, long and short path length 
transit-time variations will be treated differently. For path lengths of 
the order of hundreds of kilometers, we will assume that average-path 
rainrates of about 10 mm/hr can occur. Using a 10 mm/hr rainfall 
rate, An r is equal to 0.83 N units at 6 and 16 GHz, and 0.67 N units at 
30 GHz. The corresponding values of At' are plotted versus path 
length in Fig. 4. Comparison of Fig. 4 with Fig. 1 shows that the transit- 
time variations one expects from rain are more than an order of magni- 
tude smaller than those expected from the atmosphere itself. This im- 
plies that as far as the dynamic range of compensation equipment in 
synchronized digital systems is concerned, atmospheric variations (other 
than rain) are the determining factor for long paths. 

For short paths, rainfall rates much greater than 10 mm/hr often 
occur. Using results for New Jersey, 6 transit-time variations for short 
paths are computed using the refractive index values found in Table 
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Fig. 4 — Estimated maximum rain-caused transit-time variations for long path 
length line-of-sight communication systems. 



IV; they are plotted in Fig. 5 with P as a parameter, where P equals 
the number of minutes per year the given variation wall be exceeded. 
Variations to be exceeded 0.5 min/year and 5 min/year are shown, and 
as before, the rainfall variations are an order of magnitude smaller than 
the atmospheric variations (Fig. 2) for the corresponding path lengths. 
When viewed in terms of the nature of convective showers (showers 
exhibiting high rainfall rates), Fig. 5 aids in estimating the upper bounds 
for the time derivative of rain-caused transit-time fluctuations. Con- 
vective showers can introduce large attenuation in a path rapidly with 
onsets of the order of tens of seconds. In Fig. 5 the intersection of the 
criterion line for a 500 Mb system with the 16 GHz curve for 0.5 min/yr 
occurs at a path length of 10 km. If the onset of the rain over this path 
were to occur in 10 seconds, then i max = .01 b/s, which is a factor of 6 
smaller than the time derivative for the fronts considered previously. 



IV. CONCLUSIONS 



The maximum transit-time variations encountered in the troposphere 
are due primarily to changes in the gaseous atmosphere rather than 
rain, and may amount to 1.6 /zs for a transcontinental path; this converts 
to 800 bits in a 500 Mb system. The maximum time derivative of the 
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Table IV — Index of Refraction Variations Exceeded 
P Minutes per Year for Various Path Lengths 



(a) P = 5 min/yr. 





Path Average 








Pul h Length 


Rainfall Rate 


An r at 6 GHz 


An r at 16 GHz 


An- at 30 GHz 


(km) 


(ram/hr) 


(.V units) 


(N units) 


(A units) 


1.3 


140 


9.44 


7.42 


4.64 


2.6 


135 


9.11 


7.18 


4.51 


5.2 


110 


7.40 


5.98 


3.86 


7.8 


90 


6.14 


5.02 


3.3 


10.4 


80 


5.48 


4.54 


3.0 



(b) P = 0.5 min/yr. 





Path Average 








Path Length 


Rainfall Rate 


An- at 6 GHz 


An- at 16 GHz 


An- at 30 GHz 


(km) 


(mm/hr) 


(A* units) 


(A units) 


(N units) 


1.3 


190 


12.8 


10.1 


6.21 


2.6 


175 


11.78 


9.3 


5.72 


5.2 


150 


10.1 


7.9 


4.9 


7.8 


130 


8.78 


6.94 


4.38 


10.4 


110 


7.46 


5.98 


3.86 



Z 0.3 



S 0.2 

< 



'- O. 1 

in 

7 



30GHZ 

1 6 GHZ 

6 GHZ 



P = 0.5 




O.IxT FOR 500 Mb 
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Fig. 5 — Rain-caused transit-time variations for short path lengths. The parameter 
P equals the number of minutes/year the given variation will be exceeded. 
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transit-time variation is due to the motion of fronts, and has been derived 
herein: T max = 1.36 X 10~ 10 sec/sec which is equivalent to l/16th of a 
bit at 500 Mb. Neither At°« nor r mttX appears so large as to prohibit 
the synchronization of digital systems transmitted over line-of-sight 
propagation paths. 
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